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Abstract. We study the free expansion of a pancake-shaped Bose-condensed gas, 
which is initially trapped under harmonic confinement and containing a vortex at 
its centre. In the case of a radial expansion holding fixed the axial confinement 
we consider various models for the interactions, depending on the thickness of the 
condensate relative to the value of the scattering length. We are thus able to evaluate 
different scattering regimes ranging from quasi-three-dimensional (Q3D) to strictly 
two-dimensional (2D). We find that as the system goes from Q3D to 2D the expansion 
rate of the condensate increases whereas that of the vortex core decreases. In the Q3D 
scattering regime we also examine a fully free expansion in 3D and find oscillatory 
behaviour for the vortex core radius: an initial fast expansion of the vortex core is 
followed by a slowing down. Such a nonuniform expansion rate of the vortex core may 
be taken into account in designing new experiments. 
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1. Introduction 

Two-dimensional (2D) condensates in harmonic confinement are attracting a lot of 
attention. By varying the anisotropy parameter, measured as the ratio between the 
trap frequencies in the z- and the planar directions, flatter and flatter (pancake-shaped) 
condensates are being produced in magnetic or optical traps [TJ El E] with the ultimate 
goal of observing the special features of low dimensionality. As the gas approaches the 
2D limit, its collisional properties start to influence the boson-boson coupling parameter, 
which becomes dependent on the density of the system in contrast to the situation in 
3D where at low density the interactions are described by a constant coupling strength. 

The quantized vortex states are important in establishing the superfluid nature of 
the condensates [U 13 El- Current experiments on vortices [TUJ have utilized 
elongated geometries, but one may expect that studies on pancake-shaped condensates 
will follow. Understanding vortices in 2D is important since they are expected to play 
a role in the transition from the superfluid to the normal state [H]. Their role in 
inhomogeneous systems such as trapped gases is less well known. 

From the experimental point of view, observing a vortex in a trapped gas is difficult 
because of the smallness of the core region compared to the size of the gaseous cloud. 
Several solutions have been suggested to circumvent this problem. The size of the core 
relative to the cloud radius is larger in systems made of a smaller number of particles. 
Another possibility is to create vortices with large angular momentum as achieved in a 
recent experiment [ 12 i . The vortex core increases in a free expansion after the release 
of the trap. Indeed, it has been suggested [Hj that the vortex core may expand faster 
than the condensate cloud. Vortices were first experimentally observed by releasing the 
trap and allowing the condensate to expand [Oj. 

Several theoretical works have studied the free expansion of vortices in Bose 
condensed systems inside anisotropic trap potentials employing various 3D models for 
the coupling constant jHJ 13 EJ H2j- m a recent work [T3] we have calculated the 
equilibrium density profiles of pancake-shaped condensates using various models for 
the coupling and we have found that, when the condensate is in the 2D regime for what 
concerns its collisional properties, its density profile is very different in size from what 
is predicted by the 3D model. 

We consider in this paper pancake-shaped condensates with a vortex, employing 
various models for the coupling parameters to describe the 2D regimes of interest. We 
study both their equilibrium profile and the expansion properties after the trapping 
potential is released. By calculating the time evolution of the condensate size and of the 
vortex core radius and their expansion rates in 2D we again establish how the crossover 
regime may be identified in such experiments. We find that the condensate cloud 
expands faster in its plane as the system goes from being Q3D to 2D. The expansion rate 
of the vortex core radius, on the other hand, decreases in the dimensionality crossover. 
We also find that the expansion characteristics of the condensate cloud and of the vortex 
core are rather different when the trap potential is released only in the radial plane or 
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in both the z- and the perpendicular directions. In the latter case the velocity of the 
vortex core can be greater than that of the condensate cloud in the initial phase of the 
expansion. 

The paper is organized as follows. We first introduce our description of condensates 
with a vortex for different regimes of scattering properties in Sec. II. We then present 
our results in Sec. Ill and conclude with a brief summary in Sec. IV. 



2. The model 



We consider a dilute Bose-condensed gas at zero temperature under anisotropic harmonic 
confinement characterized by a trap frequency uj±_ in the x-y plane and an axial trap 
frequency u z = Xu±, with A ^> 1. The motion in the ^-direction is frozen and we 
may describe the ground state of the condensate by the wave function ip(r) in the 
x-y plane, to be determined from a 2D Schrodinger equation. The dynamics of the 
2D condensate which accommodates a quantized vortex state is described by a time- 
dependent nonlinear Schrodinger equation (NLSE), which reads 



di ^ ' 



h 2 h 2 K 2 l 

V 2 + + -mu 2 ± r 2 + #(</>(r,t))|^(r,t)| 2 



4>(r,t). (1) 



2m 2mr 2 2 

Here, Hk is the quantized angular momentum along the z-axis, m is the atomic mass, 
and g(if>) is the coupling strength in 2D which depends on the condensate wave function 
as discussed below. We have assumed that the order parameter entering the NLSE can 
be written as $(r, t) = i/;(r, t)e iK< ^ where is the azimuthal angle. We consider repulsive 
interactions between the atoms in the condensate, which enter the coupling strength g 
through a positive scattering length a. 

As we have discussed previously jTHj, the density profile of the condensate reflects 
the modified shape of the confining potential and the modified scattering properties. 
When the linear dimension a z = {h/muj z ) l l 2 of the condensate cloud in the ^-direction 
is much larger than the 3D scattering length (a z 3> a), the collisions take place in three 
dimensions and the coupling constant to be used in the NLSE is 

g Q3D = 2%/2T , (2) 

m a z 

including only the geometrical effects of the reduced dimensionality. 

When the anisotropy further increases and a z becomes comparable with a ( a z ~ a), 
the collisions start to be influenced by the confinement in the z-direction and the 
interaction strength assumes a form appropriate to a quasi-two-dimensional condensate, 

= 2^{h 2 /m){a/a z ) 
g Q 2D 1 + ^ nnigQ2Dn{2rrm/h 2 )a 2 )l - U 



/2na z 

Here, n(r, t) = N\if)(r, t) \ 2 is the areal density of the condensate cloud, whose dependence 
on the radial coordinate and on time is to be calculated self-consistently during the 
numerical solution of Eq. (1). This expression was originally derived by Petrov et 
al. |14| ITK] by studying the scattering amplitude in a system which is harmonically 
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confined in the ^-direction and homogeneous in the x-y plane. The same result has also 
been obtained by Lee et al. \W\ \T7\ within a many-body T— matrix approach. Note that 
the coupling strength now depends on density, as is to be expected for 2D collisions. 
Furthermore, gQ2D is given by an implicit relation which has to be solved numerically 
during the solution of the NLSE [TH] . 

Finally, the strictly 2D regime is approached as a z becomes much smaller than a 
(a z <C a) and the system is described by the coupling 

47ih 2 l 

92D — ,-; — -, — o\7 ' W 

m | m (nor) | 

as was first derived by Schick [Tj| for a homogeneous Bose gas of hard disks [2U] ■ The use 
of the coupling strength g 2 D for inhomogeneous systems, involving a dependence on the 
local density, has been proposed by Shevchenko [2IJ an d more recently by Kolomeisky 
et al. [221- Note that also g 2 D has a spatial dependence due to n but, contrary to <7q2D, 
carries no information on the confinement in the z-direction. 

We solve the time-dependent NLSE iteratively by discretization using a split-step 
Crank-Nicholson scheme j2H|. In the case of 2D calculations a simple one- dimensional 
array of grid points is sufficient to describe the spatial part of i/j(r,t) at each time 
step. The 3D calculations that we subsequently report for ip(r, z, t) as a solution of the 
standard Gross-Pitaevski equation (GPE) involve a two-dimensional grid because of the 
cylindrical symmetry of this wave function. 



3. Results and discussion 



3.1. Equilibrium profiles and radial expansion 

The numerical solution of the NLSE with a centrifugal term in equilibrium conditions 
gives the ground-state wave function of the condensate with a single vortex (k = 1). We 
first take values of the anisotropy parameter and of the scattering length as appropriate 
for 23 Na atoms in the experiment of Gorlitz et al. PP (A = 26.33 and a = 2.8 nm, so 
that a/a z = 3.8 x 1CT 3 ). We scale the radial coordinate by the harmonic oscillator 
length a± = {h/muj^) 1 ^ 2 and the radial wave function by l/a±. We take N = 5000 since 
it is easier to observe vortices for smaller numbers of particles. For these parameters, 
the system undergoes collisions in 3D but has 2D characteristics for what concerns the 
confinement effects. As in the absence of the vortex, we find that the Q3D and Q2D 
models yield wave function profiles identical to each other as shown in Fig. 1(a), whereas 
the 2D model produces a quantitatively very different result. 

We next look at the planar free expansion of the condensate with a vortex for 
the same parameters, holding fixed the confinement in the axial direction. The time- 
dependent NLSE is solved after the trap potential \mui 2 r 2 is switched off at time t = 0. 
Figure 1(b) shows the time dependence of the root-mean-square (rms) value of the radial 
coordinate r rms = (r 2 ) 1 / 2 , which describes the size of the condensate cloud, and of the 
vortex core radius R c in three different models of the coupling. As in [2], the core 
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radius is chosen to be the value of r where the condensate density reaches 1/e of its 
peak value, viz. \tp(R c )\ 2 = e _1 max{|?/>(r)| 2 }. We see from Fig. 1(b) that the strictly 2D 
model predicts a much faster expansion compared to the Q3D and Q2D models, which 
again yield identical results for the present set of parameters. The vortex core expands 
much more slowly than the condensate cloud, as indicated in the same figure. We also 
show in Fig. 1(c) the velocity of expansion of the atomic cloud and of the vortex core as 
functions of time. After an initial accelerated motion both velocities attain a constant 
value for tu± > 2. 

We increase the anisotropy parameter to A = 2 x 10 5 and take a/a z = 0.33 in 
order to investigate the regime of crossover from 3D to 2D. As shown in Fig. 2(a), the 
profiles of the wave function with a vortex all look similar in this regime of parameters. 
Consequently the time dependence of the radial coordinate and of the vortex core radius 
display similar behaviours during the expansion, as shown in Fig. 2(b). We note that 
the expansion of the condensate cloud is faster than in the Q3D regime. This is also 
indicated by the higher velocities depicted in Fig. 2(c). The velocity of the vortex core 
radius, however, decreases compared to the previous regime. The small wiggles seen in 
R c in the inset of Fig. 2(c) are due to numerical differentiation of the results shown in 
Fig. 2(b) and are not expected to be physical. 

The system enters the strictly 2D regime when we further increase the scattering 
length to a/a z = 2.68. In this case, the vortex wave function is best described by the 2D 
model using g 2 n, although the results of the Q2D model using gQ2D are very similar. As 
the initial shape of the cloud in the Q3D model extends out to a larger radial distance 
as shown in Fig. 3(a), one may again distinctly identify the truly 2D regime by the wave 
function profile. Furthermore, we observe in Fig. 3(b) that the atomic cloud expands 
slightly faster compared to the previous Q3D and crossover regimes. This is also evident 
in Fig. 3(c), where the expansion velocity of the cloud is larger. The velocity R c of the 
vortex core radius, on the other hand, is very small in the strictly 2D regime. 

The above results indicate that the initially trapped Bose condensate expands faster 
in the x-y plane as the system moves from being 3D to 2D. The opposite behaviour is 
observed for the vortex core. However, the expansion rate of the condensate cloud 
remains larger than that of the vortex in all cases that we have considered. Our 
numerical results are in good agreement with previous studies of vortices in Bose- 
condensed fluids [7J |5| for similar values of the parameters. 

3.2. Three-dimensional expansion 

We focus now on the case of a free expansion in 3D space, adopting the set of 
experimentally relevant parameters as in Fig. 1. In this case the presence of the 
confinement does not affect the binary collisions between the atoms and the system is 
well described by the 3D time-dependent GPE with coupling constant g^r, = 4irfi 2 a/m. 

We have first tested the consistence of our Q3D model by solving the GPE for 
a 3D anisotropic system to study the expansion when only the radial trap potential 
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muj 2 1 r 2 /2 is turned off {i.e. the potential muj 2 z z 2 /2 remains throughout the expansion). 
The results are identical to those obtained earlier in Figure 1 with the 2D kinematics 
and the Q3D scattering properties, showing the correctness of our physical picture. 

We have then studied the expansion properties of a 3D anisotropic condensate when 
the trap potentials in both r and z directions are released. Essentially identical results 
are obtained for two different choices of the initial configuration, i. e. starting from the 
ground state of the 3D GPE or from the ground state of the 2D GPE with Q3D coupling 
combined with a Gaussian profile in the third direction. In both cases we calculate the 
time-dependent behaviour of the rms values of the radial and axial coordinates, as well 
as that of the core radius R c . The results are shown in Fig. 4. We observe a much faster 
expansion for the z coordinate, reflecting the initially tighter axial confinement. It also 
appears that the vortex core radius undergoes a slow oscillation, as is more evident in 
Fig. 4(b) where we plot the velocities. 

The oscillatory behaviour of the vortex core can be understood as being closely 
related to the repulsive interactions, which strongly affect the density profile. As 
compared to the ground state of the noninteracting case, the repulsive interactions give 
rise to a smaller vortex core radius and to a larger rms radius. The density suddenly 
decreases when the trap is released in both r and z directions, as a result of the fast 
expansion in the z direction, and this causes the nonlinear interaction energy to rapidly 
vanish on the time scale of the radial expansion. Therefore, in the first part of the 
expansion (for tujj_ < 0.8 in our study) the compressed density at the centre of the 
cloud expands faster than the remaining part, with the result that the core expansion 
rate is larger than that of the condensate cloud. Similar results were also found by 
Dalfovo and Modugno [Zj. Then, when the density is so low that the interaction energy 
becomes negligible, the vortex cannot continue to expand faster than the cloud and the 
velocity R c even reverses sign around tu± ~ 1.5, as shown in Fig. 4(b). 

We have also studied the expansion properties of systems with N = 5 x 10 4 and 
N = 5 x 10 5 particles and the same trap frequencies and scattering length as in Fig. 1, 
finding again a much rapid expansion in the z direction and oscillatory behaviour in the 
vortex core radius. These results are not shown since they are qualitatively very similar 
to those in Fig. 4. 

4. Summary and concluding remarks 

In the perspective of experimental investigations on low-dimensional condensates, we 
have studied the free expansion of pancake-shaped Bose-Einstein condensates with a 
vortex. By choosing various values of the trap anisotropy and scattering length we 
have explored different scattering regimes, from a quasi-three-dimensional regime to a 
strictly two-dimensional one. We have described the cloud at mean field level, taking as 
coupling strengths for the boson-boson interactions the density-dependent expressions 
derived by Petrov et al. Jl] for the quasi-2D regime and by Kolomeisky et al. [22] for 
the strictly 2D regime. 
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We have considered both the case of a 2D expansion keeping the confinement along 
z fixed and the case of a fully 3D expansion. In the 2D case we have observed that with 
increasing anisotropy the expansion rate of the cloud increases while that of the vortex 
core decreases. In the 3D expansion, performed in the regime where the collisions are 
three-dimensional and the Q3D model is the most appropriate, the expansion rate of 
the vortex core exhibits an oscillatory behaviour due to the interplay between nonlinear 
interactions and anisotropic confinement. 

Whereas in this work we have restricted our analysis to the case of a condensate 
with a single vortex at T = 0, it may also be interesting to investigate the role of the 
noncondensate cloud at finite temperature, as vortices are intimately connected with 
the Kosterlitz-Thouless phase transition. To explore this area it would be necessary to 
take into account multiple vortices within the condensate cloud jHj. 
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Figure 1. Properties of a condensate with a vortex within the models described by 
9Q3D (dotted lines), <?q2D (solid lines) and g 2 D (dashed lines), for N = 5000, A = 26.33 
and a/a z = 3.8 x 10~ 3 . (a) The vortex-state equilibrium wave function, (b) The time 
dependence of the root-mean-square (rms) value of the radial coordinate r rms (upper 
curves) and that of the vortex core radius R c (lower curves), (c) The time dependence 
of the rms value v rms of the velocity of the radial coordinate (upper curves) and 
that of the velocity R c of the vortex core radius (lower curves). The dotted lines are 
superimposed to the solid lines for this set of parameters. 




Figure 2. Properties of a condensate with a vortex within the models described by 
9Q3D (dotted lines), gQ2D (solid lines) and giu (dashed lines), for N — 5000, A = 2x 10 5 
and a/a z = 0.33. (a) The vortex-state wave function at equilibrium, (b) The time 
dependence of rms radial coordinate r rms and of the vortex core radius R c (inset), 
(c) The time dependence of the rms velocity v rms of the radial coordinate and of the 
velocity R c of the vortex core radius (inset). 



Free expansion of two-dimensional condensates with a vortex 



10 




Figure 3. Properties of a condensate with a vortex within the models described by 
9Q3D (dotted lines), c/Q2D (solid lines) and giu (dashed lines), for N = 5000, A = 2x 10 5 
and a/a z = 2.68. (a) The vortex-state wave function at equilibrium, (b) The time 
dependence of the rms radial coordinate r rms and of the vortex core radius R c (inset) . 
(c) The time dependence of the rms velocity v rms of the radial coordinate and of the 
velocity R c of the vortex core radius (inset). 




Figure 4. (a) The time dependence of the rms values of the radial (solid lines) and 
^-coordinates (dotted lines), and that of the vortex core radius (dashed lines), for a 3D 
system with N = 5000, A = 26.33 and a/a z = 3.8 x 10~ 3 . (b) The time dependence of 
the velocities r rms , z rms and R c for the same system. 



